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ABSTRACT 

During a previous outburst phase, the neutron star 4U 2129+47 exhibited evidence for a spatially extended 
Accretion Disk Corona (ADC) via broad, partial X-ray eclipses occurring periodically on the 5.24 hr orbit. Since 
1983, however, 4U 2129+47 has entered a quiescent state several orders of magnitude fainter. We have performed 
a 37ksec (~2 binary orbits) Chandra observation of 4U 2129+47 to determine whether an extended coronal 
structure also exists in quiescence. Total eclipses are found, and the rapidity of the eclipse ingress and egress 
are used to place upper limits on the size of the X-ray source. The spectrum is found to be comprised of a 
soft component plus a fairly hard power law tail. The former is seen to be sinusoidally modulated over the orbital 
period in a manner consistent with neutral column variations, possibly due to the interaction of an accretion stream 
with a (small) disk about the neutron star. We fit realistic neutron star atmosphere models to the soft spectra, and 
comment on the consistency of the spectra with cooling neutron star models. It has been previously suggested that 
the 4U 2129+47 system is a hierarchical triple, with the outer body being an F star. We use differential astrometry 
to show that the X-ray point source and F star are spatially coincident to within 0.1". We further compare newly 
determined upper limits for the extrinsic neutral column to the reddening of the F star. Finally, we discuss how 
future X-ray observations can further constrain models of quiescent neutron star emission, as well as directly 
verify the triple hypothesis. 

Subject headings: accretion, accretion disks - neutron stars - stars: individual :4U 2129+47 - X-rays:stars 



1. INTRODUCTION 

Accretion Disk Corona' (ADC) sources such as 4U 2129+47 
are believed to be near edge-on accreting systems since they 
have shown binary orbital modulation via broad, partial X-ray 
eclipses. Such eclipses imply emission extended over radii 
comparable to the binary separation (McClintock et al. 1982; 
White & Holt 1982). Studies of the ADC X1822-371 offer an 
intriguing possibility for the cause of this very extended corona 
(Heinz & Nowak 2001). Orbital period derivative measure- 
ments imply that XI 822-371 may be accreting very close to 
its Eddington rate; however, only a very small fraction of the 
Eddington flux is observed. If the mass transfer from the sec- 
ondary is actually accreted by the primary, then the ADC must 
either be very optically thin (the tacit assumption made for 'on 
state' emission of 4U 2129+47 by Garcia & Grindlay 1987), or 
very optically thick. If the latter is true, the released accretion 
energy may not be able to diffuse out of the flow on a viscous 
inflow time scale, and the corona will essentially become an 
'advection dominated accretion flow' (ADAF; Esin, Narayan, 
& McClintock 1997, and references therein). In order to dis- 
sipate the accretion energy, an adiabatic outflow might ensue 
at large radius (Adiabatic Inflow/Outflow Solution', ADIOS; 
Blandford & Begelman 1999), hence creating the ADC. 

The ADC source 4U 2129+47 is strongly believed to be a 
neutron star since it has exhibited a Type I X-ray burst (Gar- 
cia & Grindlay 1987). Its classification as an ADC is due to 
the fact that prior to 1983, 4U 2129+47 was observed to have 
a lightcurve similar to that of XI 822-371. Over the 5.24 hr or- 
bit of 4U 2129+47, the eclipse width was 0.2 in phase and 



w 75% of the X-rays were occulted at the eclipse midpoint. As 
hypothesized for XI 822-371, the ADC might have been indi- 
cating the presence of an ADIOS. 

Since 1983, 4U 2129+47 has been in a quiescent state. The 
question then arises as to whether the extended coronal struc- 
ture has persisted or recurred in this state (see Garcia 1994; 
Garcia & Callanan 1999). An extended coronal structure could 
signal the presence of an ADIOS, albeit at low accretion rates 
where advection domination might ensue due to a decoupling of 
the ions and electrons. The bulk of the accretion energy would 
stored as thermal energy of the protons, which would then be 
lost as a wind rather than be dissipated on the neutron star sur- 
face (Blandford & Begelman 1999). 

Alternatively, the spectra of quiescent neutron star sources 
might be due emission from their surface as they cool after a 
prior active accretion phase, or due to a low level of residual 
accretion (van Paradijs et al. 1987; Brown, Bildsten & Rut- 
ledge 1998; Rutledge et al. 2000; Bildsten & Rutledge 2001). 
Thermal emission from the neutron star surface has been sug- 
gested as partly contributing to the putative differences between 
the fluxes of quiescent neutron stars and quiescent black hole 
candidates, with the former class being suggested as systemat- 
ically brighter (Narayan, Barret & McClintock 1997; Bildsten 
& Rutledge 2001; Garcia et al. 2001). Modern neutron star 
atmosphere models utilizing realistic opacities (Rajagopal & 
Romani 1996; Zavlin, Pavlov & Shibanov 1996) have yielded 
good fits to ROSAT spectra of quiescent neutron stars, includ- 
ing quiescent spectra of 4U 2129+47 (Rutledge et al. 2000). 
The flux and color temperature of the spectrum are determined 
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by the level of residual accretion and/or by the time aver- 
aged (over 10 4 -10 5 yr time scales) accretion rate in prior ac- 
tive phases (Brown, Bildsten & Rutledge 1998). The degree to 
which residual heat dominates over active accretion in specific 
sources is a matter of current debate, with there being a recent 
suggestion that the quiescent flux level of 4U 2129+47 might 
be too low to be consistent with the cooling model (Wijnands 
2002). 

It is also possible that for quiescent neutron stars the "pro- 
peller effect" (Campana et al. 1998, and references therein) 
halts, or at least greatly diminishes, the accretion rate onto the 
neutron star surface. The neutron star might then enter a radio 
pulsar phase, with X-ray emission arising from the interaction 
between the pulsar wind and the accretion flow from the sec- 
ondary (Campana et al. 1998; Menou et al. 1997; Campana & 
Stella 2000). The X-ray emission would be prominent at very 
large radii, comparable to the disk circularization radius. As in 
the outburst phase of 4U 2129+47, this would imply that the X- 
ray emission would be gradually eclipsed over ^ 0. 1 in orbital 
phase. 

A number of neutron star systems in quiescence have exhib- 
ited hard tails (Asai et al. 1996, 1998; Campana et al. 1998; 
Rutledge et al. 2001), which might be related to the shocked 
pulsar wind. The shocked wind emission might be in addition 
to and somewhat independent of any emission associated with 
the neutron star surface, i.e., the hard tail could vary while the 
soft excess remained steady (Campana & Stella 2000). The 
presence or absence of such a hard tail in the quiescent state of 
4U 2129+47 has been impossible to determine with ROSAT ob- 
servations (0.1-2.4keV), although those observations did sug- 
gest that the soft X-ray component was variable (Garcia 1994; 
Garcia & Callanan 1999). 

4U 2129+47 has another very interesting property which re- 
lates to estimates of its quiescent flux. After 4U 2129+47 en- 
tered quiescence, optical modulation was no longer observed 
over its 5.24 hr binary period. Instead of finding its expected 
M or K type companion, a nearby late F type star was found 
(Kaluzny 1988; Thorstensen et al. 1988; Chevalier et al. 1989; 
Garcia et al. 1989; Cowley & Schmidtke 1990), which has led 
to the hypothesis that the 4U 2129+47 system is a hierarchical 
triple. The F star would be in a few month orbit about the inner 
binary, which in turn could lead to modulations of the binary 
eccentricity on 50 yr time scales (Garcia et al. 1989). It is the 
optical measurements of this F star, which may or may not be 
related to the X-ray source, that yields the 6.3 kpc distance that 
has been used by Garcia et al. (2001) to determine the quies- 
cent flux of 4U 2129+47. We note that given a 6.3 kpc distance, 
4U 2129+47 is the third brightest quiescent neutron star in the 
Garcia et al. (2001) sample. 

The Chandra X-ray satellite, which provides excellent spa- 
tial resolution (^ 1") and good coverage in the soft X-ray band 
(0.5-10keV) can help answer some of the questions concerning 
the 4U 2129+47 system. The high spatial resolution can deter- 
mine whether or not the X-ray source is truly coincident with 
the observed F star. Measurement of the eclipse can determine 
whether an extended coronal structure still exists and place con- 
straints on the mass and radius of the secondary. The 0.5-2 keV 
spectra can yield constraints on surface emission models, while 
the 2-10keV spectra can indicate whether 4U 2129+47 exhibits 
a hard tail similar to other quiescent neutron stars. To these 
ends, on 9 December 2000, we obtained a 37 ksec (« 2 binary 
orbital periods) observation of 4U 2129+47. Details of the ob- 



servational modes and data extraction methods are presented in 
Appendix 2. 

The discussion of our results is organized as follows. We first 
outline our data extraction methodology in §2. In §3 we discuss 
the previous optical observations and the extent to which the 
X-ray source is coincident with the identified F star. We then 
present the X-ray lightcurve in §4 and discuss its implications. 
We present orbital phase resolved spectroscopy in §5. We sum- 
marize our results in §6, and suggest further observations with 
the X-ray Multi-Mirror Mission (XMM) that can resolve the 
mysteries left unanswered by these Chandra observations. So 
as not to deter the reader from the scientific results, some of 
the details of the analyses are relegated to the Appendix, in- 
cluding discussion of the statistical methods used in analyzing 
the lightcurve (Appendix A). In Appendix B, we briefly dis- 
cuss the sources that were serendipitously observed along with 
4U 2129+47. 

2. DATA ANALYSIS METHODOLOGY 

All observations of 4U 2129+47 described here were 
taken with the Advanced Camera from Imaging Spectroscopy 
(ACIS) on board the Chandra X-ray satellite. Specifically, they 
were taken with the backside illuminated S3 chip (in imaging 
mode, i.e., without the gratings inserted). To reduce the amount 
of spectral pile-up (see below), we employed a 1/4 sub-array 
on all chips, yielding a frame readout time of 1.1 sec. The total 
integration time was approximately 36.6 ksec, with roughly a 
96.5% observing efficiency (see below), yielding 35.3 ksec of 
on source time. Data extraction and response matrix generation 
were performed with the CI AO v2 . 2 analysis package and 
calibration products from CALDB v2 . 9. Sources within the 
field of view were detected using the wavdetect program, 
and source counts were extracted from within the 3cr contours 
found from this analysis. The 4U 2129+47 background region 
was extracted from an annulus centered on the source with an 
inner radius of 5" and an outer radius of 14". We include in 
our spectral and lightcurve analyses photons only in the 0.5- 
8keV range, and have rebinned the spectra to have between 
20-26 photons per spectral bin. In determining the error bars 
of the spectral parameters, we have calculated the 90% confi- 
dence limits using the methods of Lampton, Margon & Bowyer 
(1976) (i.e., Ay 2 = 2.71 for one interesting parameter). We 
have not applied any systematic errors to the spectral energy 
bins. 

Although the detected count rate from 4U 2129+47 was only 
0.05 cps, the observation did suffer from a small degree of 
pile-up (i.e., two photons being detected as a single photon with 
their summed energy). An in depth discussion of pile-up, along 
with analysis methods that may be able to account partially for 
its effects, is presented by Davis (2001). In all fits presented in 
this work, we have employed the pile-up model of Davis (2001) 
as implemented in ISIS vl.0.6 (Houck & DeNicola 2000). 

With a 1 . 1 sec frame integration time and an « 0.05 cps count 
rate, we expect 28 double photon events over the course of 
our entire observation. Roughly 95% of the detected photons 
land on the same 3x3 pixel region, which implies that 10% of 
the double photon events were registered as individual photons, 
whereas the remaining 90% of the double events were poten- 
tially registered as a single photon. The greatest uncertainty in 
this latter number is the "grade migration factor", a, i.e., the 
probability that the piled-up event is not rejected as a cosmic 
ray particle or other background event. In our use of the ISIS 
implementation of the model of Davis (2001), we have assumed 
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has 28 photons detected over the course of the observation, and 
therefore is well-localized. Applying an ss 2" translation 1 , lin- 
ear in RA and DEC, to center the X-ray image of 4U 2129+47 
on the HST image of the F star, it is found that the X-ray im- 
age of the S3-S source coincides with an optical counterpart to 
within 0.1". This is near the limit for differential astrometry 
performed with Chandra. The 4U 2129+47 X-ray source and 
the optically identified F star discussed in the references listed 
in Table Bl are therefore, to good accuracy, spatially coinci- 
dent. In §5, we comment further upon the extent to which the F 
star and the quiescent X-ray counterpart of 4U 2129+47 are at 
the same distance. 
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that there is a 50% probability (a = 0.5) that two piled-up pho- 
tons are accepted as a single "good event". 

Consistent with the above expectations, two lightcurve time 
bins show two distinguished photons. Since the detected spec- 
trum peaks at s=s 1 keV, the majority of piled-up events accepted 
as photons likely were read as having energy ^ 2keV. Even 
though we only expect < 26 piled-up events, as we discuss be- 
low there are only « 100 photons in the hard tail of our ob- 
served spectrum; therefore, pile-up does effect the fitted power 
law slope and amplitude. If a = 0.5 is in fact too low a prob- 
ability for accepting a piled-up event, then we have under sub- 
tracted the piled-up photons at 2 keV and slightly overestimated 
the value of the photon index. Likewise, if a = 0.5 is too large, 
then we have slightly underestimated the photon index, T. To 
date, there has been no independent calibration of the value of 
a (which in reality is likely to be energy and chip location de- 
pendent) in the Chandra detectors. Our estimates are that for 
the fits described below these systematic effects lead to at most 
AT « ±0.3. 

3. ASTROMETRY 

Based upon the assumption of Roche lobe overflow and a 
main sequence mass-radius relationship, the 5.24 hr orbital pe- 
riod of 4U 2129+47 suggests a secondary mass of « O.6M , 
which corresponds to an M or K dwarf. It was expected that 
optical observations of the 4U 2129+47 system in quiescence 
would reveal strong ellipsoidal variations. Instead, Thorstensen 
et al. (1988) observed a system exhibiting an F star with no de- 
tectable orbital modulations. Comparison of the position of this 
quiescent optical source with the position of the optical coun- 
terpart of 4U 2129+47 in outburst revealed the two sources to 
be coincident within 0.3"; therefore, Thorstensen et al. (1988) 
considered the possibility that the 4U 2129+47 system might 
be a triple. With higher resolution optical observations, Cheva- 
lier et al. (1989) showed that the 4U 2129+47 field was fairly 
crowded, but that there were no obvious contaminating sources. 
These authors favored the triple hypothesis. Garcia et al. (1989) 
found that the F star exhibited spectroscopic velocity modula- 
tions on the order of tens of km s" 1 on the time scale of months, 
which they also attributed to the system being a triple. 

In Table B 1 we list some of the optical identifications made 
for the quiescent counterpart to 4U 2129+47. All spectroscopic 
observations, ranging from 2200-7500 A, agree on some type 
of F star identification for this counterpart, and most observa- 
tions suggest an extinction of E(B-V) « 0.3. If the X-ray 
source is at the same distance as the F star, then it is at a distance 
of ~ 4-8 kpc. (The value of 6.3 kpc used by Narayan, Gar- 
cia, & McClintock 1997 and Garcia et al. 2001 is taken from 
Cowley & Schmidtke 1990.) With the exquisite spatial reso- 
lution of Chandra, we now can also determine to what degree 
the X-ray source is coincident with the F star. The J2000 X-ray 
source position— RA= 21 h 31 m 26.19 s , DEC= 47° 17'24.7"— 
agrees very well with the coordinates of the F star observed by 
the Hubble Space Telescope Faint Object Spectrograph (HST- 
FOS) (Deutsch et al. 1996). 

The spatial coincidence can be further tested with the HST 
images obtained (separately from the FOS spectra) by Deutsch 
et al. (1996). In addition to the image of the F star, the HST ob- 
servations also contain a likely counterpart for the source S3-S 
listed in Appendix B. The S3-S X-ray source, likely a star, 

1 As described by threads on the Chandra and HST web pages (ht tp : / / a s c . 
cross calibrate the acquired images with an accurate astrometric database, absolute astrometry errors of 0.5"-1.5" are expected for each instrument. As we are here 
only interested in differential astrometry, we have not attempted to absolutely calibrate the coordinates of each image field. 
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Fig. 1.— Folded 0.5-8 keV lightcurve of 4U 2129+47, (repeated at 
phases 1-2 for clarity). Error bars are 1-cr Poisson errors, except dur- 
ing the eclipse where the horizontal bar is the 3-cr upper limit to the 
count rate. The dashed line shows the best fitting sinusoidal modula- 
tion plus a total, rapid eclipse. 

4. LIGHTCURVE 

One of the main goals of our observation was to deter- 
mine whether or not the broad, partial dip seen in the outburst 
lightcurve of 4U 2129+47 (McClintock et al. 1982; White & 
Holt 1982) persisted into quiescence. Garcia (1994) did not 
detect any evidence of an eclipse in quiescence (see below), 
while Garcia & Callanan (1999) found only weak evidence for 
an eclipse and could not distinguish among total and partial 
eclipse models. The folded lightcurve from our observation, 
shown in Fig. 1, removes any ambiguity and shows that the 
eclipse is rapid, total, and of shorter duration than in outburst. 
The eclipse duration is 1523^0 sec (see Appendix A), i.e., 0.08 
in binary phase. This is 2.5 times narrower than in outburst. 
No photons are detected over the course of the eclipse, which 
implies that the 3a upper limit to the mean observed count rate 
during eclipse is w 0.002 cps. Given the mean count rate out of 
eclipse, the eclipsing object exhibits an optical depth of t > 3 
with 99.9% confidence, and t > 4 with 95% confidence. 

In addition to the eclipse, the lightcurve exhibits a sinusoidal 
modulation. Fitting the lightcurve outside of eclipse with a 

harvard . edu and http : / /www . st sci . edu, respectively), if one does not 
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function R((f>) = A+Bsin[27r(</>-0o)L where R is the observed 
rate and (j> is the orbital phase (<fi = is the midpoint of the 
eclipse), we find A = 0.037 ± 0.002, B = 0.012 ± 0.002, and 
0o = 0.60 ±0.03 (90% confidence levels). Note that this sinu- 
soidal modulation is similar to that exhibited by XI 822-371 in 
its bright state (Parmar et al. 2001; Heinz & Nowak 2001, and 
references therein), where it is believed that modulation is due 
to interactions of the accretion stream with the accretion disk. 
Thus the sinusoidal modulation may be indicating the presence 
of an active accretion stream and a disk about the neutron star. 

Note that Garcia (1994) in three separate, short ROSAT-HRI 
observations found evidence of variability. Due to the fact that 
these observations were spread over the span of a month, one 
could argue whether intrinsic source variability or orbital mod- 
ulation was observed. The sinusoidal modulation observed here 
is roughly consistent, in both phase and amplitude, with the X- 
ray variability described by Garcia (1994) (who likely detected 
little, if any, of the actual eclipse). Thus, whereas the quiescent 
lightcurve of 4U 2129+47 has changed from that in outburst 
(Garcia 1994), it is possible that the quiescent lightcurve has 
remained remarkably steady, at least in shape and relative am- 
plitude, during quiescence. 

With an accurate measurement of the eclipse, we can up- 
date the ephemeris presented by McClintock et al. (1982). The 
barycenter corrected Julian Date of the first eclipse (<p = 0) in 
our lightcurve is 



JD = 245 1879.5713 ±0.0002 



(1) 



which is 0.0293 days later than expected from the McClintock 
et al. (1982) ephemeris. This discrepancy could be due to either 
the uncertainty in the original orbital period determination, or 
due to orbital evolution. Adding a quadratic term to the Mc- 
Clintock et al. (1982) ephemeris, and propagating errors, we 
find the time of phase zero to be 



T Q = 2444403.7429 ±0.002 + (0.2 182579 ±8 x 10" 7 ) N 

+[(5.0±2.4)xlO" n ]A/ 2 . 



(2) 



(For these observations, we have observed orbits N = 34252 and 
34253.) Orbital evolution is indicated, but only at the 2<r level. 
Converting this into an orbital period derivative, we find 



P- 1 P = (1.6±0.8)x 10 6 years 



(3) 



Due to the possibility of winds, magnetic torques, etc., it can 
problematic to convert an observed period derivative directly 
into an implied accretion rate (see, for example, Frank, King 
& Raine 1992, for a general discussion of the role of some of 
these effects). However, if the above period derivative is real 
and solely due to accretion from the secondary to the compact 
object, the implied accretion rate averaged over the past twenty 
years is (M) = (4 ± 2) x 10 18 g s" 1 . Given that 4U 2129+47 has 
been in quiescence over most of this time, this would imply 
an extremely large accretion rate in outburst, and/or substan- 
tial mass loss without radiation, i.e., an ADIOS, in quiescence. 
Comparably large orbital period derivatives and implied mean 
accretion rates have been found for XI 822-371 (Parmar et al. 
2001; Heinz & Nowak 2001). 

We can also determine an ephemeris for the average period. 
Specifically, we have 

r = 2451879.5713±0.0002+(0.21825962±6xl0- 8 )JV , (4) 

where the period at the time of these observations might be as 
much as 10" 6 days longer than the above average period. 



0.012 



0.010 



0.008 



•§ 0.006 
o 



0.004 



0.002 



0.000 



\ 
A 

. \ \ 









- 


; V- 


\ 




- 




\ \ 








\ \ 








\ \ 








\ \ 








\ \ 








\\ 








\\ 








V 


















x^ 













50 100 150 

Transition Width (s) 



200 



250 



Fig. 2. — Bayesian probability distribution for the width of the 
eclipse ingress (dashed line) and egress (solid line), assuming a linear 
transition between the eclipsed and uneclipsed states. 

The fact that the eclipses observed with Chandra are total 
and only 0.08 wide in phase already indicates that the X-ray 
emission is emanating from a region nearly two times smaller 
than in outburst. The rapidity with which eclipse ingress and 
egress occurs, however, can place even tighter constraints. As- 
suming that the eclipse is total and that eclipsing medium has a 
sharp edge, then the time it takes for the count rate to go from 
the mean count rate down to zero, or the reverse, is indicative 
of the size of the X-ray source. As described in Appendix A, 
we have performed a statistical analysis of the eclipse ingress 
and egress, and show the results in Fig. 2. We have averaged 
the two ingresses with each other, and the two egresses with 
each other, but we have not averaged the ingresses with the 
egresses. At the 95% confidence level, the average ingress is 
more rapid than 170 sec, while the average egress is more rapid 
than 217 sec. This implies that the X-ray source diameter is 

0.07 a, where a is the binary separation. For a total system 
mass of 2M Q , this implies an X-ray diameter of < 9 x 10 9 cm 
(approximately 7 earth diameters). This is approximately a fac- 
tor of 15 smaller than the implied X-ray source diameter for 
4U 2129+47 in outburst (McClintock et al. 1982; White & Holt 
1982). 

An important caveat for the above estimates is that they per- 
tain to the 0.5-8 keV lightcurve only, and are dominated by 
emission from ^ 2keV. Thus, we are unable to distinguish be- 
tween different source sizes above and below 2keV, as have 
been hypothesized (Campana et al. 1998; Campana & Stella 
2000). We also note that for the most part, each individual 
ingress and egress looked statistically identical, with the ex- 
ception of the second egress. There was week evidence for 
a prolonged period of decreased emission towards the begin- 
ning of the second egress. This was most readily noted by the 
fact that the \a limit for the second eclipse egress was 260 sec, 
as opposed to 130-170 sec for the other ingresses and egress. 
Accretion "dips", due to "blobs" of enhanced neutral hydro- 
gen column are in fact common near zero phase in Cyg X-l, 
for example (see Balucihska-Church et al. 1997 and references 
therein). Perhaps a similar occurrence is being hinted at here, 



Total Eclipses of 4U 2129+47 



5 



although the statistics are not strong enough to say with any 
degree of certainty whether or not this is the case. 

The lightcurve not only places constraints upon the X-ray 
source diameter, but the measured eclipse duration also places 
constraints upon the parameters of the secondary, specifically 
its mass and radius 2 Here we calculate the mass ratio of the 
secondary to the primary, q, as a function of inclination angle, 
;. Assuming that the secondary completely fills its Roche lobe, 
which we approximate as spherical, then the ratio of the eclipse 
duration to the orbital period is given by 



(5) 



where a is the orbital separation of the binary and R2 is the 
Roche lobe radius of the secondary. Utilizing the approxima- 
tion of Paczyhski (1971), this can be rewritten as q rts qo(l —qi) 
where 




qo 



2 -1 3 /2 

+ cos z ;J 



^[2cos 2 /-(? 2 ] [e 2 + cos 2 i 



1/2 



(6) 



The other interesting aspect to note about the observed 
lightcurve is that the position of the (out of eclipse) lightcurve 
minimum occurs at a completely different phase than in out- 
burst. In the lightcurve presented by McClintock et al. (1982), 
the (weak) non-eclipse minimum occurred at phase ~ 0.8. That 
is, the dip leads the eclipse by w 0.2 in phase, similar to what 
is found for XI 822-371 (Parmar et al. 2001; Heinz & Nowak 
2001). In quiescence, the non-eclipse minimum occurs at phase 
~ 0.4, i.e., it leads the eclipse by 0.6 in phase. In a model where 
the minimum corresponds to the location of where the accretion 
stream first interacts with the accretion disk, this shift can be 
explained if the disk has substantially shrunk in size. 

As a rough illustration of this, consider a crude model where 
the incoming matter has the specific angular momentum solely 
due to orbital motion at the inner Lagrange point, and other- 
wise is only falling in the gravitational potential of the compact 
object (i.e., we are assuming q <C 1). Assuming an elliptic or- 
bit intersecting the disk, it is relatively straightforward to show 
that the radius of the accretion disk, R D , relative to the binary 
separation is given by 



Assuming specific values for the primary mass (1.3 and 
2.2 Mq), in Fig. 3 we plot the secondary mass, M2, and Roche 
lob radius, R2, as a function of orbital inclination. 
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where x = (a - R2) /a, R2 is the Roche lobe radius of the sec- 
ondary, and (f> is the phase (in radians) by which the dip leads 
the eclipse. (The fact that here <j> > n indicates that the approxi- 
mations are not fully valid; however, they are sufficient to illus- 
trate the general point.) Approximating the size of the Roche 
lobe as a function of mass ratio as above, one finds, for a broad 
range of q, that R D « 0.4-0.5 a for <j> = 0.2. This is consistent 
with independent estimates of the disk radius of XI 822-371 
(Heinz & Nowak 2001). Taking cj) = 0.5 and q = 0.2, one finds 
the substantially smaller value of Rd « 0.2a. More accurate 
estimates of the disk radius would rely upon obtaining better 
estimates of the binary inclination and mass ratio, q. 

5. SPECTRUM 

In Fig. 4 we present the 0.5-8 keV spectrum of 4U 2129+47, 
averaged over all phases except the eclipse. Note that in this 
figure we have not corrected for photon pile-up, although we 
do account for this effect in our fits. Overall, the 0.5-2 keV 
spectrum is very similar to the ROSAT Position Sensitive Pro- 
portional Counter (PSPC) spectrum described by Garcia & 
Callanan (1999). The Chandra spectrum, however, clearly 
shows the presence of a hard tail, even after accounting for pile- 
up, that lies beyond the upper energy cutoff of ROSAT. 

As a first description of the Chandra spectrum, we fit a model 
consisting of a blackbody plus power-law that are absorbed by a 
neutral column (utilizing the absorption model of Wilms, Allen, 
& McCray 2000). Results of this fit are presented in Table B2 
(model A), where we have parameterized the blackbody nor- 
malization by an effective radius, R, equivalent to assuming 
an emitting area of nR 2 and a source distance of 6kpc. The 
power-law flux is for the photon index being oc E~ r , and the 
power law normalization is photons keV" 1 cm" 2 s" 1 at 1 keV. 
(Note the large error bars on this normalization are more in- 
dicative of the uncertainties in the power law slope, as opposed 

2 Previous estimates of these parameters have been presented by Home, Verbunt & Schneider (1986) utilizing measurements taken during outburst of a sinusoidally 
modulated H/3 absorption line and a double peaked He II 4686 emission line. The two measurements together implied a primary and secondary mass of 0.6 ± 0.2Mq 
and ().4±O.2M0, respectively. These estimates, however, are hampered by uncertainties in modeling the outburst emission of the secondary and disk, as well as by 
any possible contamination from the F star. 
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Fig. 3. — Constraints on the mass (solid lines) and radius (dashed 
lines) of the secondary, as a function of the system inclination. We 
have assumed a primary mass of 1.3 Mq (lower lines) and 2.2 M 
(upper lines). 

If one takes as a straw-man model that the primary is a neu- 
tron star with mass 1.3-2.2M© (as evidenced by the Type I 
burst observed by Garcia & Grindlay 1987) and that the Roche 
lobe radius of the secondary is greater than or equal to its main 
sequence radius (as has been measured for the secondary in 
numerous accreting systems), then the eclipse measurements 
alone imply that the system inclination i > 75° and that for the 
secondary mass 0.25 M < M 2 < 0.6 M Q (i.e., 0.2 < q < 0.5). 
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to the flux uncertainties at energies > 2keV.) The neutral col- 
umn, blackbody temperature, and implied flux we obtain are 
consistent with the fits to PSPC observations described by Gar- 
cia & Callanan (1999), which were also averaged over phase. 
The flux and spectrum of 4U 2129+47 is therefore consistent 
with having remained constant since 1994. 3 




1 8 
Energy (keV) 



Fig. 4.— Measured 0.5-8 keV spectrum of 4U 2129+47, excluding 
the eclipses, fitted with an absorbed blackbody plus power law. Lower 
panel shows the residuals from this fit. 

The fitted radius of < 4 km (90% confidence levels) for these 
observations is significantly smaller than the w 10 km expected 
for a neutron star; however, as discussed by Rutledge et al. 
(2000), this might be attributed to the fact that a simple black- 
body is not the correct model for a realistic neutron star atmo- 
sphere. Numerous theoretical models exist that more carefully 
treat the radiative transfer in such an atmosphere. To describe 
the Chandra spectra of 4U 2129+47, we here use an ISIS 
(Houck & DeNicola 2000) implementation of the model of Za- 
vlin, Pavlov & Shibanov (1996). Fit parameters for this model 
include neutron star temperature, mass, radius, and distance. 
For all fits described here, we have fixed the neutron star mass 
to 1.4M Q . As presented in Table B2, we have performed fits 
with a fixed distance of 6kpc (model C), and fixed radii of 
10km (model E) and 5 km (model G). The parameter space for 
these fits contain a number of local \ 2 minima, and the "best fit" 
parameters and their associated error bars describe the statistics 
within the parameter regions of these local minima. 

We see that fixing the distance at 6 kpc leads to good fits with 
a slightly larger neutron star radius of > 5 km and a temperature 
very similar to that of the simple blackbody fit. Postulating a 
larger neutron star radius of 10 km yields a lower temperature 
of « 0. 1 keV — consistent with that found by Rutledge et al. 



(2000) using these same models to describe the PSPC data; 
however, the inferred distance is then rather small. The 90% 
confidence upper limit of < 4 kpc is less than most of the dis- 
tance estimates based upon optical studies of the F star. (Note 
that the fit with the fixed 10 km radius also yielded the largest 
fitted value of the neutral column.) Assuming instead a neu- 
tron star radius of 5 km yields distance estimates consistent with 
those of the F star. 
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Fig. 5. — Phase resolved spectra of 4U 2129+47 (upper panel). 
The filled circles represent the 2-8 keV spectrum, measured from all 
phases excluding the eclipses. Diamonds and asterisks represent the 
spectrum integrated over phases 0.65-0.96 and 0.2-0.5, respectively. 
Dashed line is the best fit absorbed blackbody plus power law model 
where all parameters, except neutral hydrogen column, were con- 
strained to be identical as a function of phase. Lower panel shows 
the contribution to the x 2 of the residuals from this fit. 

It had already been noted by Garcia & Callanan (1999) that 
the neutral column obtained from the orbital phase-averaged 
ROSAT PSPC spectrum was too large to be consistent with 
the reddening of the F star. Using our typical fitted column 
of 4 x 10 21 cm -2 implies E(B-V)« 0.75 (Predehl & Schmitt 
1995), whereas the E(B-V)= 0.34 measured by Deutsch et al. 
(1996) implies a column of 1.9 x 10 21 cm" 2 . Phase resolved 
spectroscopy offers a possible solution to this discrepancy. As 
shown in Fig. 1, there is clear phase-dependent variability. Ex- 
tracting spectra at orbital phases 0.2-0.5 and 0.65-0.958, the 
variation appears to be consistent with changes in the soft spec- 
trum only. (The statistics of the hard tail above 2keV become 
poor in each phase bin, but are consistent with a constant ampli- 
tude and slope.) We therefore have chosen to perform a series 
of fits where below 2keV we consider the above two orbital 
phases, while above 2keV we use the spectrum from the to- 



3 Earlier ROSAT High Resolution Imager (HRI) observations (Garcia 1994) lacked any detailed spectral information. However, assuming the same model as for 
the PSPC fits, Garcia & Callanan (1999) argued that during the 1992 HRI observations, 4U 2129+47 was approximately a factor of two more luminous. Rutledge 
et al. (2000), using neutron star atmosphere models, also argued that 4U 2129+47 was likely more luminous during the HRI observations. 
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tal observation, excluding the eclipse. We constrain all model 
parameters, including normalizations, to be the same for each 
of these spectra, except that we allow the neutral column to be 
different in each orbital phase. (For the hard tail, we constrain 
the neutral column to be the same as for orbital phases 0.65- 
0.958.) All three spectra are fitted simultaneously, and the re- 
sults are presented in Table B2 and Fig. 5. As for our orbital 
phase-averaged spectra, we present a simple blackbody model 
(model B), and neutron star atmosphere models with fixed dis- 
tances (model D) or radii (models F, H). 

Overall these fits are fairly good, and we see that the fitted 
temperatures, radii, distances, and power law parameters are 
very comparable to those obtained for the phase-averaged fits. 
The neutral column is seen to vary by nearly a factor of two 
between the two phases, with the lower column value being 
a factor of only 10-70% larger than would be implied by the 
F star reddening. The averages of the two column values are 
approximately the same as the columns determined from the 
phase-averaged fits. The spectral variations are therefore con- 
sistent with being solely due to variations of the intrinsic col- 
umn, with the total column having roughly equal contributions 
from intrinsic and extrinsic columns. 

As regards the power law components in these fits, all fits 
imply that the power law comprises <; 10% of the (unabsorbed) 
0.5-8 keV flux. The error bars are large; however, the index 
of the power law is consistent with being fairly flat. We note 
that, if anything, changing our assumptions regarding pile-up is 
most likely to further flatten the power law (see §2). In terms 
of amplitude, slope, and extent, such a power law is difficult to 
describe solely using a neutron star atmosphere model. 

6. DISCUSSION 

Our Chandra observations of 4U 2129+47 have revealed that 
the previously observed extended coronal structure has col- 
lapsed by at least a factor of 15 in size. The soft X-ray spectrum 
is now consistent with thermal emission from the surface of a 
neutron star with radius on the order of 5 km. To what extent 
this emission is indicative of residual thermal emission or ac- 
tive, albeit low-level, accretion, remains unclear. The fact that 
the spectrum observed by Chandra is consistent in amplitude 
and shape with that observed by ROSAT in 1994 can easily 
be accommodated for in the cooling neutron star model, as the 
cooling time scale is expected to be very long. The prior ob- 
served soft X-ray variability (Garcia 1994) can be accounted 
for via orbital phase-dependent obscuration effects, rather than 
intrinsic accretion variability. 

On the other hand, questions concerning the cooling model 
remain. First, assuming a distance of 6kpc, consistent with 
various estimates of the F star distance, the fitted neutron star 
radius tends to be be fairly low, on the order of 5 km. Although 
larger than the inferred blackbody radius, this is smaller than 
expectations from numerous neutron star models. Forcing a 
large neutron star radius, 10 km, implies a distance substan- 
tially less than most estimates of the F star distance, which of 
course might still be a chance overlap. However, it is seen that 
each of the various fits, including the simple blackbody, are sta- 
tistically nearly identical to one another. Second, the variabil- 
ity discussed by Garcia (1994) is consistent, in amplitude and 
phase, with that observed here, which indicates that roughly 2/3 
of the integration time from the 1992 HRI observation occurred 
over the phases of lowest emission. Thus the estimate of an ap- 
proximate factor of two flux difference between 1992 and 1994 
might be an underestimate, although this is hampered by the 



poor spectral resolution of the HRI. Rutledge et al. (2000) ar- 
gue that this change might be related to solely the outer crust 
cooling. Additionally, this change in flux might indicate that as 
of 1992, 4U 2129+47 had not yet truly entered into a quiescent 
phase. 

Based upon the work of Brown, Bildsten & Rutledge (1998) 
and the observations discussed by Garcia (1994), Garcia & 
Callanan (1999), and Rutledge et al. (2000), Wijnands (2002) 
has suggested that the quiescent phase of 4U 2129+47 might be 
too faint to be consistent with the neutron star cooling model. 
Specifically, one expects a quiescent flux of 



where (F ) is the average flux in outburst, t is the average time 
the source spends in outburst, and t q is the average time spent in 
quiescence. Assuming a typical accretion efficiency of e ~ 0.1, 
and an average accretion rate of (M) ~4x 10 18 g s" 1 , as sug- 
gested by the possible period derivative found in §4, we then 
find that 




To accommodate the above result, either the assumed accretion 
rate is atypically large (i.e., not indicative of the long term aver- 
age in outburst), the durations of the quiescent periods are very 
long, there is enhanced cooling of the neutron star during quies- 
cence (see Wijnands 2002, Wijnands et al. 2002, and references 
therein), or the efficiency is very low. 

If one instead assumes an average outburst flux of (F q ) ~2x 
10~ 9 erg cm" 2 s" 1 (i.e., 5 times the typical 2-10keV flux seen in 
outburst; Garcia & Grindlay 1987 and references therein), then 
one obtains t q /t„ 20. However, this estimate assumes that 
one has observed nearly all of the outburst flux in the system, 
contrary to the more usual assumptions for an ADC of either 
a very optically thin or very optically thick coronae seen near 
edge on (Parmar et al. 2001 ; Heinz & Nowak 2001). In order to 
retain the typical ADC assumptions, one is again drawn toward 
hypotheses of the previous outburst episode being atypically 
bright, long quiescent periods, enhanced neutron star cooling, 
or very low radiative efficiencies. Given that any matter that 
reaches the neutron star surface will contribute towards heating 
the neutron star, the latter hypothesis implies that inefficiency 
during outburst is achieved via an outflow, e.g., an ADIOS. 

Questions remain about the origin of the sinusoidal variabil- 
ity and the power law tail. The former is similar to what has 
been observed in XI 822-371, where the variability has been 
associated with the interaction of the accretion stream with the 
outer edge of an accretion disk (Parmar et al. 2001; Heinz & 
Nowak 2001). A similar model for 4U 2129+47 could argue for 
active accretion. Observations of XI 822-371 imply that a to- 
tally opaque disk rim is raised that completely blocks the central 
neutron star, while a more vertically extended disk atmosphere 
(with Nh 10 22 cm" 2 ) attenuates the extended corona (Heinz & 
Nowak 2001). Here the evidence suggests that the neutron star 
in 4U 2129+47 is directly observable, while the combination of 
inclination angle and outer disk atmosphere height serves only 
to attenuate its soft flux with a column of Nh ~ 2 x 10 21 cm" 2 . 
The presence of the hard tail might also argue for active accre- 
tion. The limited Chandra band pass shows that at least > 10% 
of the total, unabsorbed 0.5-8 keV luminosity emanates from 
this component. Thus we might be seeing the accretion energy 
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first being funneled through a corona, and then actively being 
dissipated as thermal energy on the neutron star surface. 

One alternative that preserves the cooling model would be for 
the accretion stream to be halted by an outgoing pulsar wind 
(Campana et al. 1998; Menou et al. 1997; Campana & Stella 
2000). The hard tail would represent shocked emission from 
the interaction region of the wind and accretion stream. The si- 
nusoidal modulation could still be due a neutral hydrogen col- 
umn being raised by the interaction of the accretion stream with 
a disk rim. Such a disk need not be actively accreting, however, 
as the pulsar wind and magnetic field could halt any inward flow 
of disk material. We note that the disk radius inferred from the 
location of the light curve minimum is comparable to expecta- 
tions from the pulsar wind model; however, the soft flux radius 
determined from the rapidity of the eclipse is somewhat smaller. 

The Chandra observations have lent further credence to the 
hypothesis that the 4U 2129+47 system is a hierarchical triple. 
The X-ray source position is now seen to be even more closely 
aligned with optical position of the F star. Furthermore, the 
minimum fitted neutral column (i.e., the flux from the maxi- 
mum emission phase of the lightcurve) is more consistent with 
the reddening of the F star than were previous measurements 
obtained from phase-averaged spectra. 

Prospects are good for refining the answers to some of the 
remaining questions about 4U 2129+47. Further observations 
with either Chandra or XMM-Newton, aside from improving 
the statistics of the model fits, might reveal variability. Variabil- 
ity in the soft spectra would argue against the cooling model be 
the sole source of quiescent emission. XMM-Newton observa- 
tions can have direct bearing upon the pulsar wind models for 
the hard tail. Aside from extending detection of this tail to be- 
yond 8 keV (which would rule out most atmosphere models for 



the hard tail), the methods of §4 then could be applied sepa- 
rately to the hard tail. Measurement of three or more eclipses 
would allow one to place limits of « 0.1 a on the diameter of 
the hard component emission region. Although large, this is 
a smaller diameter than suggested by many of the pulsar wind 
models. 

Additionally, each XMM-Newton observation can determine 
the eclipse centroid to approximately ±8 sec. If the F star is 
truly in a month long orbit about the binary, then we expect 
±80 sec variations in the eclipse ephemeris. Thus, multiple 
XMM-Newton observations spaced out over a month long time 
scale could directly determine the orbital parameters of a third 
body. In the absence of the detection of such orbital varia- 
tions, either Chandra or XMM-Newton observations could de- 
termine whether or not the orbital period evolution hinted at in 
the analysis of §4 actually occurred. This in turn, as discussed 
above, could have implications for the expected level of quies- 
cent emission in the cooling neutron star model. 
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ECLIPSE ANALYSIS 

In what follows, we explicitly ignore the effects of detector "deadtime". Each Chandra frame for this observation represented an 
integration time of approximately 1.1s with a "frame transfer time" of 40 ms between frames 4 . The effective deadtime therefore was 
on the order of 3.5%, which is negligible for the purposes of this analysis. We shall employ a Bayesian approach. 

The following discussion pertains to eclipse egress; however, it is readily generalized to ingress as well. We label our At = 1.1s 
time bins with indices — > k, and assume that at the beginning of the time bin indexed by 0, the X-ray source is fully eclipsed, and 
by the beginning of the time bin indexed by k it is fully unobscured. We assume that the mean count rate measured at the detector 
increases linearly from to r (although other assumptions are readily accommodated). The mean count rate in time bin i is therefore 
r, = r(i+ l/2)/k, for i = — > k—l. Following the usual derivation of Poisson statistics, we assume that each bin is short enough 
such that the probability of measuring one photon event in a given bin is P(l) = r,Af, the probability of measuring no events is 
P(Q) = 1 — P(l), and that P(> 1) = 0. The probability of the data, D, given this model is therefore 



P(D|r,fc) = []r,Af[](l- Af) 



(Al) 



where the product over ; is for those bins where a photon was detected, and the product over j is for those bins where no photons 
were detected. We divide the above by the probability assuming a steady, uneclipsed count rate of r, which, given r <C 1, yields 



P(D\r,k) cx J] (7) l[[l + (r-rj)At] = J] 



1+ 1/2 



n 



1 + rAf 



k-j-l/2 



(A2) 



Further approximating the product over empty time bins as the exponential of a sum, taking a uniform prior for k (the eclipse 
sharpness) and a gaussian prior for the detected count rate, and employing Bayes theorem, we obtain the probability distribution for 
the model parameters as 



P(r,k\D) oc exp [-(r-rf /2a 2 r ] x JJ 



i+l/2 



exp 



k-j-l/2 



(A3) 



4 See the Chandra Proposer's Observatory Guide, http://cxc/harvard.edu. 
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In practice, based upon the fits to the folded lightcurve, we have used r = 0.046 cps and 0.042 cps for the mean ingress and egress 
count rates, respectively, and we have taken o> = 0.0035 cps for both. In the above we have assumed a fixed start point for the egress 
commencement. We allow for uncertainty in this start point by assuming that its prior probability is uniformly distributed between 
and 500 s before the first detected photon. Eq. (A3) then must be slightly modified by including in the sum within the exponential all 
the time bins from 500 s before the first detected photon to the time bin represented by index 0. 

The above was used to calculate probability distributions separately for the two measured eclipse ingresses and the two measured 
eclipse egresses. The probability distributions were then multiplied together for each pair to give the ingress and egress distributions 
for the folded lightcurve. Marginalizing over the eclipse sharpness yielded the eclipse widths and midpoints, while marginalizing 
over the eclipse start points yielded the distributions for the eclipse sharpness (see §4). 

SERENDIPITOUSLY OBSERVED SOURCES 

A number of other sources were serendipitously observed within a few arcminute radius of 4U 2129+47. These sources, along 
with their coordinates and 0.3-8 keV total counts (i.e., integrated over the entire 36.6 ksec observation) are presented in Table B3. 
For the most part, these sources can be identified with stellar objects in the Digital Sky Survey 5 . Several sources, specifically S3-/3, 
S2-a, I2-a, and I3-a might be faint active galactic nuclei. S3-/3 contains sufficient counts to fit a crude spectrum, and it is found 
to be consistent with an absorbed power law, with A^h = 3 ± 3 x 10 21 cm -2 and V = 1 .9 ± 0.6, 90% confidence limits. Its 0.5-8 keV 
(absorbed) flux corresponds to 4 x 10" 14 erg cm" 2 s" 1 . There is also weak evidence that its count rate increases by a factor of two on 
time scales of « 25 ksec. 

The source S3-<5 was useful for performing differential astrometry, as it likely was detected also with the HST observations of 
Deutsch et al. (1996). For purposes of differential astrometry (as opposed to Table B3), the source position was determined with the 
program wavdetect applied to the 0.5-2 keV image. 

Note that given the broad energy range (0.3-8 keV), and large confidence contours (3-<r detection region from wavdetect), the 
total counts in Table B3 might contain photons unrelated to the source, and therefore should be regarded as upper limits. This is 
especially true as the angular distance from the target point, S3-a, increases. 
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Table B 1 



Optical identifications for the quiescent counterpart to the 4U 2129+47 system. 



Wavelength 


E(B-V) 


Dist. 


Spec. Type 


Reference 


4500-7500 A 


0.3 


> 1 kpc 




Thorstensen et al. 1988 


4800-6600 A 


0.5 


4kpc 


F5-8 V 


Chevalier et al. 1989 


4000-4800 A 






F7IV-V 


Garcia et al. 1989 


3700-5600 A 


0.3 


6-8 kpc 


F8IV 


Cowley & Schmidtke 1990 


2200-3300 A 


0.34 




F8 IV 


Deutsch et al. 1996 



Table B 2 

Fits to the X-ray Spectra of 4U 2129+47 (90% Confidence Level error bars). 



Fit 


N H 


kT 


R or D a 


A r 


r 


0.5-2keV b 


0.5-2 keV c 2- 


8keV 


X 2 /DoF 




(10 21 cm- 2 ) 


(keV) 


(km or kpc) 


(10- 6 ) 




(10- 


' 13 erg cm" 2 s" 1 ) 






A 
B 


o 1+0.9 
J-1 -Q.7 

2 +0 ' 9 

z - u -0.8 
^• J -1.0 


21 +0 02 

u - zl -0.03 

22 +0 03 

u - zz -0.02 


2 5 +L5 

z - 1 -0.7 


6.4 +20 2 2 
2.8!' 4 


0.6!};i 


1.0 
1.3 
0.8 


2.5 
2.5 


0.5 
0.5 


38.1/33 
37.3/31 


C 


o 4+O.3 


2S+ 00 
u - z;, -0.06 


dg q+0.9 


? 1 +0 - 5 

z - 1 -1.8 


0.4 +0 f 


1.0 


2.9 


0.5 


38.2/33 


D 


9 7 +0.4 
z - '-Q.7 

5 0+ 06 


0.24^ 


ls.o +2 - 4 


1 5 +5 - 2 

1 - J -1.3 




1.3 
0.8 


3.1 


0.5 


38.0/31 


E 


4 2 +L2 

^• z -0.9 


0.08to1 


1 9 +21 


4 6 +201 
^•°-3.9 




1.0 


4.0 


0.5 


37.6/33 


F 


J - z -l.l 
5 4 +L4 


1 0+ 002 


2 7 +3 - 3 

z -'-1.6 


9 9+8.I 
z ' z -2.0 




1.3 
0.8 


4.0 


0.5 


37.9/31 


G 


T, Q+1.0 
J '^-0.9 


22 4<) 05 

u - zz -0.05 


4.2! 3 « 


3.7 + ' 5 3 8 


0.8+J 2 


1.0 


3.6 


0.5 


37.4/33 


H 


2 - 8 -o:8 


?s +0 - 05 

u - z;5 -0.05 


«r 0+4^2 
J-O-3.0 




0.3t& 


1.3 


3.6 


0.5 


38.0/31 




5 1 +L3 


0.8 









"Neutron star radius (km) assuming a distance of 6 kpc (A-D), or neutron star distance (kpc) assuming a radius of 
10km (E-F) or 5 km (G-H). 

b Absorbed Flux. 

c Unabsorbed Flux. 

d Lower limit set by table boundaries of neutron star atmosphere model. 



Table B3 
Sources in the Chandra field. 



Source 11 


RA 


DEC 


Counts b 




(J2000) 


(0.3-8 keV) 


S3-a 


21 h 31 m 26.19 s 


47° 17' 24.7" 


1301 


S3-/3 


21 h 31 m 24.57 s 


47° 17' 12.7" 


205 


S3-7 


21 h 31 m 30.29 s 


47° 16' 16.5" 


36 


S3-6 


21 h 31 m 27.86 s 


47° 17' 32.5" 


28 


S3-e 


21 h 31 m 29.21 s 


47° 16' 36.9" 


25 


S2-a 


21 h 3 l m 34.64 s 


47° 15' 3.6" 


88 


S2-/3 


21 h 31 m 33.18 s 


47°15'4.3" 


9 


YI-ol 


21 h 30 m 55.13 s 


47° 09' 49.6" 


115 


I3-a 


21 h 30 m 27.83 s 


47°15'23.5" 


80 



a Chip ID-source identifier (brightest to faintest). 

b Counts from within the 3-cr detection region as deter- 
mined via wavdetect. 



